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Legged locomotion on flowing ground (e.g. granular media) is un- 
like locomotion on hard ground because feet experience both solid- 
and fluid-like forces during surface penetration. Recent bio-inspired 
legged robots display speed relative to body size on hard ground 
comparable to high performing organisms like cockroaches but suf- 
fer significant performance loss on flowing materials like sand. In 
laboratory experiments we study the performance (speed) of a small 
(2.3 kg) six-legged robot, SandBot, as it runs on a bed of granu- 
lar media (1 mm poppy seeds). For an alternating tripod gait on 
the granular bed, standard gait control parameters achieve speeds 
at best two orders of magnitude smaller than the 2 body lengths/s 
(si 60 cm/s) for motion on hard ground. However, empirical adjust- 
ment of these control parameters away from the hard ground set- 
tings, restores good performance, yielding top speeds of 30 cm/s. 
Robot speed depends sensitively on the packing fraction </> and the 
limb frequency uj, and a dramatic transition from rotary walking to 
slow swimming occurs when 4> becomes small enough and/or lu large 
enough. We propose a kinematic model of the rotary walking mode 
based on generic features of penetration and slip of a curved limb in 
granular media. The model captures the dependence of robot speed 
on limb frequency and the transition between walking and swim- 
ming modes but highlights the need for a deeper understanding of 
the physics of granular media. 
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Compared to agile terrestrial organisms, most man-made 
vehicles possess limited mobility on complex terrain [1] 
and are easily thwarted by materials like rubble and sand. In- 
creased locomotive performance of engineered platforms de- 
mands better understanding of interaction with complex en- 
vironments. At the same time, there is increasing evidence 
that small legged machines can have greater maneuverability 
than large wheeled vehicles in many natural environments ^ . 
However, while wheeled and treaded locomotion on sand has 
been well studied by pioneers like Bekker [3], study of the 
interaction of animals or legged devices with complex media 
like sand is in its infancy in part because the physics of 
penetration and drag in granular media is largely unexplored 
for realistic conditions. Nearly all previous experiments and 
models of terrestrial locomotion were developed for running 
and walking on rigid, flat, non-slip substrates in which the 
possibility of substrate flow was ignored O [6] [71 [8] |9] . 

Rainforest, grassland, polar tundra, mountains, and desert 
are examples of complex Earth terrains with flowing sub- 
strates that challenge locomotors; the limited experience of 
the Mars Rovers supports the presumption that extraterres- 
trial landscapes will be even more daunting. Deserts, common 
in nature and occupying about 10 percent of land surface on 
Earth ilOJ , consist largely of granular media, a representative 
complex substrate. Granular materials, defined as collections 
of discrete particles, can exhibit solid-like [11] behavior below 



a critical yield stress [15], while fluid-like |13) . gas- like [14) . 
and even glass-like [15] behaviors are possible during flow. 
Yet, compared to other complex materials like debris, mud, or 
snow, granular materials are simple enough that fundamen- 
tal understanding of the collective physics can be achieved 
through interplay of experiment and theory. Unlike more het- 
erogeneous real-world environments, granular media can be 
precisely controlled using laboratory scale devices [151 116] to 
create states of varying material strength that mimic differ- 
ent deformable flowing materials produced during locomotion 
on complex terrains. Here we systematically explore the per- 
formance of a small legged device, SandBot, on granular me- 
dia prepared in different packing states with volume fraction 
ranges typical of desert sand \^\. Despite SandBot 's (and 
its predecessor RHex [18]) ability to move nimbly and rapidly 
over a wide range of natural terrain, we find that on granu- 
lar media the locomotion is remarkably sensitive to substrate 
preparation and gait characteristics, which points to both the 
need for a more sophisticated understanding of the physics 
of motion within granular media and the possibility of bet- 
ter robotic design and control paradigms for locomotion on 
complex terrains. 

Results and Discussion 

The robot we study, SandBot (Fig. la), is the smallest 
(mass 2.3 kg) in a successful series of biologically inspired ]19] 
hexapedal robots, the RHex class [18.. RHex incorporates 
the pogo-stick like dynamics observed in a diversity of biolog- 
ical organisms running on hard ground 20]. This dynamics, 
called the Spring Loaded Inverted Pendulum (SLIP) template, 
is hypothesized to confer passive self-stabilization properties 
to both biological and robotic locomotors TT . RHex was 
the first legged machine to achieve autonomous locomotion 
at speeds exceeding one body length/s [18] and it and its 
"descendants" such as EduBot/SandBot, Whegs [22] and iS- 
prawl [23] are still the leaders in legged mobility (roughly, 
speed and efficacy) on general terrain. In fact, prior to the 
recent development of the much larger BigDog [23] platform 
(1 m long, 75 kg), RHex remained the only class of legged ma- 
chine with documented ability to navigate on complex, nat- 
ural, outdoor terrain of any kind and has been used as the 
standard platform in comparisons with commercial wheeled 
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and tracked vehicles like Packbot [251 . 

SandBot moves using an alternating tripod gait in which 
two sets of three approximately c-shaped legs rotate syn- 
chronously and TV out of phase. A clock signal (Fig. Ic), 
defined by three gait parameters (see Materials and Meth- 
ods), prescribes the angular trajectory of each tripod. The 
c-legs distribute contact [26] over their surface and allow the 
robot to move effectively on a variety of terrain. On rigid, no- 
slip ground SandBot's limb trajectories are tuned to create a 
bouncing locomotion J_8j that generates speeds up to 2 body 
lengths/s (« 60 cm/s). We tested this clock signal on granular 
media, but found that the robot instead of bouncing adopts 
a swimming gait in which the legs always slip backward rel- 
ative to the stationary grain bed and for which performance 
is reduced by a factor of 30 to ~ 2 cm/s (see Supporting In- 
formation Movie SI). We surmised that this was due to an 
interval of double stance (both tripods in simultaneous con- 
tact with the ground), which is useful on hard ground during 
bouncing gaits but apparently causes tripod interference on 
granular media. Changing the gait parameters to remove the 
double stance allowed SandBot to move (see Movie S2) in the 
granular media at speeds up to 1 body length/s (« 30 cm/s) 
in a rotary walking gait that resembles the pendular gait of 
the robot on hard ground [27| but with important kinematic 
differences (discussed below). No amount of gait parameter 
adjustment produced rapid bouncing locomotion on granular 
media. We hypothesize that the ~ 50% decrease in top speed 
relative to hard ground is associated with the inability of the 
robot to undergo the arial phases associated with the bounc- 
ing gait [28j . 

In the desert, animals and man-made devices can en- 
counter granular media ranging in volume fraction from (j) = 
0.55 to (j) = 0.64 [13, and some desert adapted animals (like 
lizards) can traverse a range of granular media with little loss 
in performance [IS]. To test the robot performance on con- 
trolled volume fraction granular media, we employ a 2.5 m 
long fluidized bed trackway (Fig. lb) , which allows the 
flow of air through a bed of granular media, in this case 
~ 1 mm poppy seeds. With initial fluidization followed by 
repeated pulses of air |30] , we prepare controlled volume frac- 
tion states with different penetration properties [31) . In this 
study, we test the performance (forward speed Vx) of SandBot 
with varied leg-shaft angular frequency (lu) for volume fraction 
{(j)) states ranging from loosely to closely packed (0 = 0.580 
to (j> = 0.633). We chose forward speed as a metric of perfor- 
mance since it could be readily measured by video imaging. 
We hypothesized that limb frequency would be important to 
robot locomotion since the substrate yield strength increases 
with volume fraction and the yield stress x robot limb area 
divided by the robot mass x velocity is proportional to the 
maximum limb frequency for efficient locomotion. 

We find that robot performance (speed) is remarkably sen- 
sitive to (see Movie S3). For example, at tj = 16 rad/s the 
robot speed Vx{t) shows a change in average speed Vx of nearly 
a factor of five as (j) changes by just 5 % (Figs. ld,e). For the 
closely packed state {(j> — 0.633), Vx ~ 20 cm/s with 5 cm/s 
oscillations during each tripod rotation, whereas for a more 
loosely packed state (0 = 0.600), Vx ~ 2 cm/s with 1 cm/s 
oscillations in velocity. 

This sensitivity to volume fraction is shown in the aver- 
age robot speed vs. volume fraction (Fig. le). For fixed u), Vx 



is effectively constant for <j) above a critical volume fraction 
<j)c{Lo), but is close to zero for below 4>c(!-^)- For fixed u, 
<j)c{i-o) separates volume fraction into two regimes: the "walk- 
ing" regime (i^ > <^c, >> 0) and the "swimming" regime 
{(j) < (f>c{uj), Vx ~2 cm/s). See Movies S4 and S5 for examples 
of rotary walking and swimming modes. 

The rotary walking mode is dominant at low u and high 
<j). In this mode, a tripod limb penetrates down and backward 
into the ground until the granular yield stress exceeds the limb 
transmitted inertial, gravitational, and frictional stresses at a 
depth d(u,(f>). At this point, rather than rolling forward like 
a wheel, the c-leg abruptly stops translating relative to the 
grains and begins slipping tangentially in the circular depres- 
sion surrounding it; at the same time, the center of rotation 
moves from the axle to the now stationary center of curvature 
(see Fig. 3a). The simultaneous halt in both vertical and 
horizontal leg motion is apparently due to the large reduction 
in friction forces which occurs when the weight of the robot is 
supported by the limbs rather than the underside of the body. 
The ensuing rotary motion propels the axle and consequently 
the rest of the robot body along a circular trajectory in the 
X — z plane with speed Ru, where R = 3.55 cm is the c-leg 
radius. The forward body motion ends when, depending on cj) 
and uj, either the second tripod begins to lift the robot or the 
underside of the robot contacts the ground. 

With increased uj limbs penetrate further as the requisite 
force to rapidly accelerate the robot body to the finite limb 
speed {Ruo) increases. As the penetration depth approaches 
its maximum 2R — h, where h = 2.5 cm is the height of the 
axle above the flat underside of the robot, the walking step 
size goes to zero since there is no longer a point in the cycle 
where the limb ceases its motion relative to the grain bed. 
Any subsequent forward motion is due solely to thrust forces 
generated by the swimming-like relative translational motion 
of the limb though the grains. Note that (/>c(tj) increases with 
UJ, and that the transition from rotary walking to swimming 
is sharper in Vx for higher uj and smoother for lower uj. The 
much slower swimming mode occurs for all volume fractions 
for UJ > 28 rad/s. 

Plotting the average robot speed as a function of limb fre- 
quency (Fig. 2a) shows how the robot suffers performance 
loss as its legs rotate more rapidly. For flxed <j>, Vx increases 
sub-linearly with a; to a maximal speed Vx at a critical limb 
frequency ujc, above which Vx quickly decreases to « 2 cm/s 
(swimming) 32 . The robot suffers partial performance loss 
(mixed walking and swimming) below ujc and near total per- 
formance loss (pure swimming) above ujc- Performance loss 
for <j> ^ 0.6 is more sudden {Auj ~ 1 rad/s) compared to 
performance loss for (j) ^ 0.6. Both ujc and Vx display tran- 
sitions at « 0.6 (Figs. 2b,c). The transition a.t (j) ~ 0.6 
for the rapidly running robot is noteworthy since it has been 
observed that granular media undergo a transition in quasi- 
static penetration properties at ~ 0.6 ^31) . 

Starting with the observed kinematics of rotary walking 
with circular slipping, we constructed a straightforward two- 
parameter model that captures the essential elements deter- 
mining granular locomotion for our legged device and agrees 
well with the data (dashed lines in Fig. 2a). The model, which 
incorporates simplifled kinematics and granular penetration 
forces while still agreeing well with a more realistic treatment 
(for a more detailed discussion of the model, see Materials 
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and Methods) , indicates that reduction of step length through 
increased penetration depth is the cause of the sub-linear in- 
crease in Vx with uj and the rapid loss of performance above 
Wc- The model assumes that the two tripods act indepen- 
dently, that the motion of each tripod can be understood by 
examining the motion of a single c-leg supporting a mass m 
equal to one-third of SandBot's total mass, and that the un- 
derside of the robot rests on the surface at the beginning of 
limb/ground contact. 

Using the geometry of rotary walking (see Fig. 3a), the 
walking step size per c-leg rotation is s = 2^y — (d + h — RY, 
where d is the maximum depth of the lowest point on the leg. 
After the robot has advanced a distance s, the body again 
contacts the ground and the c-leg moves upward. Since dur- 
ing each clock signal period there are two leg rotations (one 
for each tripod), the average horizontal velocity is 2s x limb 
frequency or = ^ . The maximum limb penetration depth 
d is thus the key model component as it controls the step 
length (see Fig. 3b) and consequently the speed. Maximum 
limb penetration depth is determined by balancing the vertical 
acceleration of the robot center of mass ma with the sum of 
the vertical granular penetration force kz |33) and the gravita- 
tional force mg, where g is the acceleration due to gravity, and 
fc((/)) is a constant characterizing the penetration resistance of 
the granular material of volume fraction 0. 

At small UJ, ma ~ = J]] = mg — kd so d = mg/k, 
which is the minimum penetration depth. For finite ui the 
penetration depth is greater since an additional force must be 
supplied by the ground to accelerate the robot body to the 
leg speed Ru when the c-leg stops translating in the material. 
Taking a — Av/At, with Aw — Ruo — and At the charac- 
teristic elastic response time of the limb and grain bed, gives 
the acceleration magnitude a = Ruo/ At. The direction of the 
acceleration depends on the position of the c-leg. To keep the 
model simple we approximate the vertical component of the 
acceleration with its magnitude. Equating the vertical forces 
with mass x acceleration (see Fig. 3c), —m^ = mg — kd, 
gives c-leg penetration d = -^("sf d) with average horizon- 

tal velocity = ^S^^l - + *) + ^ - l] Fits 

to this model are indicated by dashed lines in Fig. 2b. The 
expression captures the sub-linear increase in with lu at 
fixed k{(j}), the increase in speed at fixed u as the material 
strengthens (increasing k with increasing 0), and transition 
to zero rotary walking velocity when uu is sufficiently large. 

The expression for Vx is determined by the two fit param- 
eters k and At. The parameter k characterizing the pene- 
tration resistance increases monotonically with <j) from 170 to 
220 N/m and varies rapidly below (j> ~ 0.6 and less rapidly 
above. Its average value of ~ 200 N/m corresponds to a shear 
stress per unit depth of a ~ 470 kN/m^ (using leg area=wR 
where w is the leg width) which is in good agreement with 
penetration experiments we performed on poppy seeds that 
yield a = 300 and 480 kN/m^ for </> = 0.580 and 0.622 respec- 
tively and is comparable to previous measurements of slow 
penetration into glass beads |34) where a ~ 250 kN/m''. In 
contrast. At varies little with <j) and has an average value of 
0.4 s compared to the robot's measured hard ground oscilla- 
tion period of 0.2 s when supported on a single tripod. The 
difi'erences in At can be understood as follows. In our model 



we assume the two tripods do not simultaneously contact the 
ground; however, in soft ground this is not the case, which 
consequently reduces the effective step size per period from 
2s to a lesser value. The fit value of At is sensitive to this 
variation; reducing the step size (and thus the speed) in the 
experimental data by just 13 % decreases At to 0.2 s while k 
is increased by less than 10 %. 

Our model indicates that for deep penetration the walking 
step length is sensitive to penetration depth {e.g. Fig. 3b). 
As the walking step size goes to zero with increasing uj or de- 
creasing (j), the fraction of the ground contact time that the leg 
slips through the grains (swimming) goes to one. Swimming 
in granular media differs from swimming in simple fluids as 
the friction dominated thrust and drag forces are largely rate 
independent at slower speeds |33l 135) . When thrust exceeds 
drag and using constant acceleration kinematics, the robot 
advances a distance proportional to the net force divided by 
uj^ per leg rotation, and, consequently, speed is proportional 
to uj~^ . This explains the weak dependence of Vx on uj in the 
swimming mode. The increase in robot speed with decreas- 
ing u} is bounded by the condition that the robot speed in a 
reference frame at rest with respect to the ground cannot ex- 
ceed the horizontal leg speed in a reference frame at rest with 
respect to the robot's center of mass. This condition ensures 
the existence of and eventual transition to a walking mode as 
LU is decreased. 

The transition from walking to swimming appears gradual 
for (f) ^0.6 since the penetration depth increases slowly with u 
at small uj (Ruj/At <^ g) and the uj~'^ contribution to the per 
cycle displacement from swimming is relatively large (see e.g. 
the data at = 12 rad/s in Fig. 3b). However, for (j) ^ 0.6, 
the transition is abrupt. This sharp transition occurs because 
the step size is reduced sufficiently that the legs encounter ma- 
terial disturbed by the previous step; we hypothesize that the 
disturbed material has lower (j) and k. At higher 0, the volume 
fraction of the disturbed ground is significantly less than the 
bulk which increases penetration and consequently greatly re- 
duces s. This is not the case for the transition from walking to 
swimming at lower (and low uj) where the volume fraction 
of the disturbed material is largely unchanged relative to its 
initial value. For the robot to avoid disturbed ground it must 
advance a distance R on each step, i.e. s > R, or in terms 
of the penetration depth, d < (-^ -I- 1)_R — h = 5.0 cm (green 
dashed lines in Fig. 3b, c). The disturbed ground hypothesis is 
supported by calculations of the step size derived from the av- 
erage velocity 2s = 2ttvx/uj which show a critical step size near 
s/R = 1 at the walking/swimming transition (Fig. 3d). The 
somewhat smaller value of s/R ~ 0.9 evident in the figure can 
be understood by recognizing that for s slightly smaller than 
R the majority of the c-leg still encounters undisturbed mate- 
rial. Signatures of the walking/swimming transition are also 
evident in lateral views of the robot kinematics (see Movies 
S3-S5). 

At higher uj in the swimming mode, limbs moves with 
sufficient speed to fling material out of their path and form 
a depression which reduces thrust because the limbs are not 
as deeply immersed on subsequent passes through the mate- 
rial. However, as limb speed increases further, thrust forces 
becomes rate dependent and increase because the inertia im- 
parted to the displaced grains is proportional to uj'^. Between 
strokes, the excavated depression refills at a rate dependent 
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on the difference between tiie focai surface angle and the an- 
gle of repose [36) and the depression size. Investigating the 
competition between these different processes at high lo and 
their consequences for locomotion could be relevant to under- 
standing how to avoid becoming stranded or to free a stranded 
device. 

Conclusions 

Our studies are the first to systematically investigate the per- 
formance of a legged robot on granular media, varying both 
properties of the medium (volume fraction) and properties of 
the robot (limb frequency). Our experiments reveal how pre- 
carious it can be to move on such complex material: changes 
in (j) of less than one percent result in either rapid motion or 
failure to move, and slight kinematic changes have a similar 
effect. A kinematic model captures the speed dependence of 
SandBot on granular material as a function of 4> a^nd uj. The 
model reveals that the sublinear dependence of speed on uj 
and the rapid failure for sufficiently small 4> and/or large lj are 
consequences of increasing limb penetration with decreasing <j) 
and/or increasing u, and changes to local (j) due to penetration 
and removal of limbs. While detailed studies of impact and 
penetration of simple rigid objects exist [371 133) , further ad- 
vances in performance (including increases in efficiency) and 
design of limb geometry will require a more detailed under- 
standing of the physics associated with penetration, drag, and 
crater formation and collapse, especially their dependence on 
4>. Better understanding of this physics can guide development 
of theory of interaction with complex media advanced enough 
to predict limb design [38] and control [39l strategies, similar 
to the well-developed models of aerial and aquatic craft. Anal- 
ysis of physical models such as SandBot can also inform loco- 
motion biology in understanding how animals appear to move 
effortlessly across a diversity of complex substrates [401 126). 
Such devices will begin to have capabilities comparable to or- 
ganisms; these capabilities could be used for more efficient 
and capable exploration of challenging terrestrial {e.g. rub- 
ble and disasters sites) and extra-terrestrial {e.g. the Moon 
and Mars) environments. 

Materials and Methods 

Limb kinematics SandBot's six motors are controlled by a clock sig- 
nal to follow the same prescribed kinematic path during each rota- 
tion and, as shown in previous work on RHex, changes in these kine- 
matics have substantial effects on robot locomotor performance 41 . 
The controlling clock signal consists of a fast phase and a slow phase 
with respective angular frequencies. The fast phase corresponds 
to the swing phase, and the slow phase corresponds to the stance 
phase. A set of three gait parameters uniquely determines the clock 
signal configuration: 9s, the angular span of the slow phase; 9o, the 
leg-shaft angle of the center of the slow phase; and dc, the duty 
cycle of the slow phase. Specifying the cycle average limb angular 
frequency u) fully determines the limb motion. 

In pilot experiments we tested two sets of clock signals: a hard 
ground clock signal (HGCS) with (e^ = 0.85 rad, do = 0.13 rad, 
dc = 0.56) which generates a fast bouncing gait (60 cm/s) on hard 
ground 1181 but very slow (~ 1 cm/s) motion on granular media, 
and a soft ground clock signal (SGCS) with (9s = 1.1 rad, do = -0.5 
rad, dc = 0.45) which produces unstable motion on hard ground 
but regular motion on granular media. These experiments showed 
that the locomotor capacity of SandBot is sensitive to the clock sig- 
nal. Careful observation of limb kinematics revealed that the hard 
ground clock signal fails on granular media because of the simul- 
taneous stance phase of two tripods. In this study, we use SGCS 
and explore robot performance as a function of limb frequency and 



substrate volume fraction. 

Integrated motor encoders record the position and current (and 
thus torque) of SandBot's motors vs. time. Comparison of the 
measured and prescribed angular trajectories for both sets of tim- 
ing parameters show a high degree of fidelity with an error of a 
few percent. Therefore, SandBot's change in performance between 
HGCS and SGCS timing comes from the physics of the substrate 
interaction. 

Trackway volume fraction control To systematically test Sand- 
Bot's performance vs. substrate volume fraction, we employ a 2.5 m 
long, 0.5 m wide fiuidized bed trackway with a porous plastic 
(Porex) flow distributor (thickness 0.64 cm, average pore size 90 
/im). Four 300 L/min leaf blowers (Toro) provide the requisite air 
flow. Poppy seeds are chosen as the granular media because they are 
similar in size to natural sand 42 and are of low enough density to 
be fluidized. The air flow across the fluidized bed is measured with 
an anemometer (Omega Engineering FMA-900-V) and is uniform 
to within 10 percent. 

A computer controlled fluidization protocol sets the volume 
fraction and thus the mechanical properties of the granular me- 
dia. A continuous air flow initially fluidizes the granular media in 
the bubbling regime. The flow is slowly turned off leaving the gran- 
ular media in a loosely packed state ((/> = 0.580). Short air pulses 
(ON/OFF time = 0.1/1 s) pack the material 1301 . Increasing the 
number of pulses increases up to a maximum oi 4> = 0.633. Vol- 
ume fraction is calculated by dividing the total grain mass by the 
bed volume and the intrinsic poppy seed density. The mass is mea- 
sured with a precision scale (Setra). The density of the granular 
media is measured by means of displacement in water. In exper- 
iment, since the horizontal area of the fluidized bed trackway is 
flxed, volume fraction is set by controlling the height of the granu- 
lar media {e.g. volume = area X height). 

Kinematics measurements To characterize SandBot's motion, 
we record simultaneous dorsal and lateral views with synchronized 
high speed video cameras (AOS Switzerland) at 100 frames/s. The 
center of mass (dorsal landmark) and the axles of the right-side 
front and rear motors (lateral landmarks) are marked with reflec- 
tive material (WhiteOut). A rail-pulley system allows the robot's 
power and communication cables to follow the robot as it moves to 
minimize the drag from the cables. For each trial, we prepare the 
trackway with the desired volume fraction and place the robot on 
the prepared granular media at the far end of the trackway with 
both tripods in the same standing position. An LED on the robot 
synchronizes the video and robot motor encoder data. After each 
trial, MATLAB (The MathWorks) is used to obtain landmark co- 
ordinates from the video frames and calculate Vx{t). Three trials 
were run for each combination of (0, uj) that was tested. 

Detailed discussion of rotary walking locomotion model The 
model presented in the main body of the manuscript simplifles the 
underlying physics while capturing the essential features determin- 
ing robot speed. Here we describe a more complete model (which 
lacks a simple expression for v^) and compare its predictions to 
those of the simple model. The exact expression for the vertical 
acceleration component of the body when the limbs gain purchase 

is maz = ma sin 9 = ma^J^^ ^^^^ — ^ ''j^^ ^ instead of the approx- 
imation maz = ma used in the simple model. Using the exact 
expression, the vertical granular force necessary for walking still 
has the same peak value of m{a + g) but decreases to mg when the 
limb is at its lowest point. 

The second approximation we employed in the simple model is 
that the grain force on the leg is kz. This expression is only strictly 
valid for a flat bottomed vertically penetrating intruder 33 . Since 
the leg is a circular arc, the leg-grain contact area and the vertical 
component of the grain force are functions of limb depth and leg- 
shaft angle. Generalizing kz to a local isotropic yield stress given 
by az |12| . the vertical force on a small segment of the limb Rdip 
in length at depth z is dF^ = wazRdip cosip, where w is the limb 
width and ip the angular position of the segment with respect to 
a vertical line passing through the axle. The total vertical compo- 
nent of the force acting on the leg is then Rwa z cos ■(/>iii/>. 
Substituting z{ip) = R{cosip — 1) + d and integrating gives = 

Rwa ^^{ip + costp sinip) + {d — R) sintf^ , where ipmax = 
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cos~^(l — d/Tt) and Vmin = ipmax when the leg tip is above the 

center of the c-leg and cos^^ (^i^ ^ ^ when it is below the 
center of the c-leg. is the angular extent of the limb beyond tt 
{e.g Al; = for a semi-circular limb). 

Fig. 4a shows that the full model using realistic parameters 
shares the same essential physics as the simple model. For a given 
material strength (blue curves), the penetration depth increases 
with increasing cj (intersection of blue and red curves) until ei- 
ther the step size is reduced below the critical value (vertical green 
dashed line) or the granular force no longer provides the required 



force to support the robot's mass and accelerate it. Fig. 4b presents 
fits to the experimental data of the average speed Vx vs. lo for the 
full and simple models for Vx < v* at each cp. The fits and fit pa- 
rameters for the simple (At = 0.4 s, a = 470 kN/m^) and full 
(At = 0.2 s,a = 330 kN/m^) models are in good agreement when 
the step size is less than the critical value s = R. 
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List of Figures 

1 Locomotion of a legged robot on granular media is sensitive to substrate packing and limb frequency, (a) The six- 
legged robot, SandBot, moves with an alternating tripod gait (alternate triplets of limbs rotate tt out of phase); 
arrows indicate members of one tripod, (b) Pulses of air through the bottom of the fluidized bed trackway 
control the initial volume fraction of the granular substrate; air is turned off before the robot begins to move, 
(c) Tripod leg-shaft angle 8 vs. time is controlled to follow a prescribed trajectory with two phases: a slow stance 
phase and a fast swing phase. Overlapping trajectories from trials with — 0.633 (red) and — 0.600 (blue) at 
w = 16 rad/s demonstrate that the controller maintains the desired kinematics independent of material state, (d) 
Identical tripod trajectories produce different motion for (j) — 0.600 and — 0.633. (e) For given limb frequency 
{u} = 4, 8, 12, 16, 20, 24, and 30 rad/s) the robot speed is remarkably sensitive to (p. Red and blue circles show 

the corresponding states in (c) and (d) |^ 

2 Average robot speed vs. limb frequency, (a) For a given volume fraction (p, Vx increases sub-linearly with oj to a 
maximal average speed v'^ at a critical limb frequency Uc above which the robot swims {vx ~ 2 cm/s). The solid 
lines and symbols are for (p = 0.580, 0.590, 0.600, 0.611, 0.616, 0.622, and 0.633. The dashed lines are fits from 
a simplified model discussed in the text. (b)(c) The dependence of luc and v* on (p shows transitions at (/> « 0.6 
(dashed lines) |^ 

3 (a) Schematic of a single robot leg during a step in granular media. After reaching penetration depth d, the leg 
rotates about its center and propels the robot forward a step length s. The solid shape denotes the initial stage 
of the rotational motion and the dashed shape indicates when the limb begins to withdraw from the material 
(end of forward body motion), (b) Step length vs. penetration depth (blue) with critical step size (green dashed 
horizontal) and critical penetration depth (green dashed vertical) indicating where the robot begins to encounter 
ground disturbed by the previous step, (c) Granular penetration force for k = 1.75, 2.00, 2.25, 2.50, 2.75 x 10^ N/m 
(blue) and force required to initiate rotary walking for uu — 0,8, 16, 24, 32 rad/s (red) vs. penetration depth using 
simplified walking model with At = 0.2 s. The penetration depth at constant (p is determined by the intersections 
of the corresponding blue line with the red lines. Beyond the critical depth (green dashed line) limbs encounter 
disturbed material and move to lower blue lines, (d) Step length as a function of uj derived from 2s — 2ttv/uj 
reveals the condition for the onset of swimming as s/R « 1 [l^ 

4 (a) Non-dimensionalized granular force (blue curves) for a x 10~^ — 2.50, 2.75, 3.00, 3.25, and 3.50 g cm~^s~^ 
and the required force to initiate rotary walking, a^/g + 1 (red curves) for lj = 0, 5, 10, 15, 20, 25 and 30 rad/s 
for the full model as a function of limb penetration depth with a 225° c-leg arc angle and At = 0.15 s. The 
intersection of the red and blue curves determines the penetration depth of the limb and consequently the step 
size. At constant material strength (blue) d increases with increasing u), while at constant u; increasing material 
strength reduces d. The vertical green dashed line indicates the critical penetration depth beyond which the leg 
encounters material disturbed by the previous step, (b) Comparison of Vx vs. lo for simple (red dotted curve) 
and full (blue solid curve) models. Models are fit to the measured robot speed (symbols) for Vx < wj. The green 
dashed line indicates Vx = Rui/n or equivalently s — R. In both figures h = 2.5 cm, R — 3.55 cm, w = 1.2 cm 

and m = 767 g [ill 
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Fig. 1. Locomotion of a legged robot on granular media is sensitive to substrate packing and limb frequency, (a) The six-legged robot, SandBot, moves with an alternating 
tripod gait (alternate triplets of limbs rotate 7V out of phase); arrows indicate members of one tripod, (b) Pulses of air through the bottom of the fluidized bed trackway 
control the initial volume fraction (p of the granular substrate; air is turned off before the robot begins to move, (c) Tripod leg-shaft angle 6 vs. time is controlled to follow 
a prescribed trajectory with two phases: a slow stance phase and a fast swing phase. Overlapping trajectories from trials with = 0.633 (red) and (p = 0.600 (blue) at 
a; = 16 rad/s demonstrate that the controller maintains the desired kinematics independent of material state, (d) Identical tripod trajectories produce different motion for (p 
— 0.600 and (f> — 0.633. (e) For given limb frequency (a; = 4, 8, 12, 16, 20, 24, and 30 rad/s) the robot speed is remarkably sensitive to (p. Red and blue circles show the 
corresponding states in (c) and (d). 
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Fig. 2. Average robot speed vs. limb frequency, (a) For a given volume fraction (p, Vx increases sub-linearly with a; to a maximal average speed at a critical limb 
frequency lUq above which the robot swims {Vx ^ 2 cm/s). The solid lines and symbols are for (p = 0.580, 0.590, 0.600, 0.611, 0.616, 0.622, and 0.633. The dashed lines 
are fits from a simplified model discussed in the text. (b)(c) The dependence of lOc and on (f) shows transitions at (p ^ 0.6 (dashed lines). 
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Fig. 3. (a) Schematic of a single robot leg during a step in granular media. After reaching penetration depth d, the leg rotates about its center and propels the robot forward 
a step length s. The solid shape denotes the initial stage of the rotational motion and the dashed shape indicates when the limb begins to withdraw from the material (end 
of forward body motion), (b) Step length vs. penetration depth (blue) with critical step size (green dashed horizontal) and critical penetration depth (green dashed vertical) 
indicating where the robot begins to encounter ground disturbed by the previous step, (c) Granular penetration force for A: = 1.75, 2.00, 2.25, 2.50, 2.75 X 10^ N/m (blue) 
and force required to initiate rotary walking for a; = 0, 8, 16, 24, 32 rad/s (red) vs. penetration depth using simplified walking model with At = 0.2 s. The penetration 
depth at constant (j) is determined by the intersections of the corresponding blue line with the red lines. Beyond the critical depth (green dashed line) limbs encounter disturbed 
material and move to lower blue lines, (d) Step length as a function of a; derived from 2s = lifv j lo reveals the condition for the onset of swimming as s/i? 1. 
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Fig. 4. (a) Non-dimensionalized granular force (blue curves) for a X lO"** = 2.50, 2.75, 3.00, 3.25, and 3.50 g cm^^s""^ and the required force to initiate rotary 
walking, a^/y + 1 (red curves) for lo = Q, b, 10, 15, 20, 25 and 30 rad/s for the full model as a function of limb penetration depth with a 225° c-leg arc angle and 
Ai = 0.15 s. The intersection of the red and blue curves determines the penetration depth of the limb and consequently the step size. At constant material strength (blue) 
d increases with increasing UJ, while at constant UJ increasing material strength reduces d. The vertical green dashed line indicates the critical penetration depth beyond which 
the leg encounters material disturbed by the previous step, (b) Comparison of Vx vs. iO for simple (red dotted curve) and full (blue solid curve) models. Models are fit to 
the measured robot speed (symbols) for Vx ^ '^x- '^^^ green dashed line indicates IF* = Rio/n or equivalently s = R. In both figures h = 2.5 cm, R = 3.55 cm, 
w = 1.2 cm and m = 767 g. 



Footline Author 



PNAS I Issue Date | Volume | Issue Number | 11 



